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Abstract 

Effective critical exponents for the Xcp 4 scalar field theory are calculated as 
a function of the renormalization group block size and inverse critical 
temperature f3 c . Exact renormalization group equations are presented up to 
first order in the derivative expansion and numerical solutions are obtained 
with and without polynomial expansion of the blocked potential. For a finite 
temperature system in d dimensions, it is shown that (5 C = (3 c k determines 
whether the d-dimensional (/3 C <C 1) or {d + l)-dimensional {(3 C ^> 1) fixed 
point governs the phase transition. The validity of a polynomial expansion of 
the blocked potential near criticality is also addressed. 
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The calculation of critical exponents using renormalization group (RG) methods has a 
long standing history. For the scalar A0 4 theory in four dimensions it has been shown that 
only a trivial Gaussian fixed point exists, whereas for 2 < d < 4 a strong-coupling Wilson- 
Fisher fixed point appears [[!]]. Studies of the transition from four- to three-dimensional 
systems typically rely on the e expansion [fj]; however, through years of labor it has been 
shown that this expansion is asymptotic, requiring complicated resummations or outright 
truncation of the series. In the last few years a new tool has emerged for studying phase 
transitions in arbitrary dimension: the numerical and/or analytical solution of exact renor- 
malization group equations When solving RG equations one conventionally Taylor 
expands the blocked potential and solves for the flow of the Taylor coefficients; however, we 
find evidence that Taylor expansion of the blocked potential near criticality is not justified. 
In this Letter, we will present results obtained solving finite temperature RG equations 
numerically without polynomial expansion. 

For generality we first consider a finite temperature system defined on the manifold 
S 1 x R d with the radius of S 1 equal to the inverse temperature (3. As T — > 0, the radius 
of S l goes to infinity making the manifold equivalent to R d+1 , while in the opposite limit, 
T — > oo, the manifold becomes equivalent to R d . However, as we will show below, the critical 
behavior of a finite temperature system is actually parameterized by (3 C = f3 c k Q , where k Q is 
the renormalization blocking scale. The critical behavior of the system will be characteristic 
of (d + l)-dimensions when f3 c ^> 1 and <i-dimensions when (3 C <C 1, so that in the true 
thermodynamic limit, k o ^-0, the system will exhibit (^-dimensional critical behavior for all 
critical temperatures T c > 0. 

The bare action for a scalar field theory is given by 

S[*] = [drjd d x {y (<9.0) 2 + ^(V0) 2 + V{</>)} , (1) 

where we have set % = 1, and Z T and Z s are the tree-level temporal and spatial wavefunc- 
tion renormalization constants, respectively. Note that two wavefunction renormalization 
constants are required at finite temperature due to the breaking of Lorentz invariance. 
The RG evolution of the finite temperature blocked action 5/3^ [$] can be obtained from 

e -s,, k m= /L>[0]n5(0 M (x,r)-$(x,r))e-^, (2) 

J X 

where ^^(x, r) is the blocked quantum field, which in general can be expressed in momen- 
tum space as 

fe ,n(p, uj n ) = p[%(p, uJ n ) <KP, U n ) , (3) 

with (p, w„) being the blocking function for the system. We choose the blocking function 
to be 

p < j?l(p,u n ) = 5 n=0 e(k-\p\), (4) 

so that only the modes with |p| > k and n = are included in the blocked field with the 
other modes being integrated over in the path integral. As for the spatial component of the 
blocking function, a sharp cutoff cleanly separates the high- and low-momentum modes and 
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FIG. 1. Temperature dependence of the thermal mass parameter /jZ, coupling constant A^, and 
the expectation value (4>)p near T = T c . The T dependence of Xp in the symmetric phase is also 
illustrated. 



does not introduce any dependence on the shape of the blocking function, as is the case for 
a smooth cutoff. 

In order to derive the coupled RG equations for the blocked potential and wavefunction 
renormalization we evaluate (f2[) at one- loop level and expand the background field as $(x) = 
$ + $(x). Derivative terms are handled by treating the fluctuation field, $(x), and the 
momentum, p, as operators obeying the commutation relation pi, $(x) = i<9j$(x) ||. 

We note here that since $(x) is independent of r we can neglect the RG flow of Z T and 
make the notational identification Z = Z s throughout the rest of the Letter. In terms of the 
bare potential V{Q) and bare wavefunction renormalization Z the one-loop contributions to 
the blocked potential and wavefunction renormalization are 




(5) 

where m p> $ = (Z p 2 + V" ^($)) 1 ^ 2 , and A is the ultraviolet cutoff and a tilde indicates a 
perturbative quantity. 

Due to their perturbative nature, these equations will fail to capture the true physics 
near criticality. This is because as the infrared (IR) cutoff, k, is lowered, the loop expansion 
of the path integral will no longer be expanding about the true minimum of the blocked 
action but instead about the minimum of the bare action. As an improvement, we can take 
advantage of the arbitrariness of the IR cutoff and divide the integration into a large number 
of thin shells. For instance, we could integrate eq. (|5|) from the UV cutoff, A, down to A — Ak 
and use the potential obtained, uj^\_^ k , as the initial condition for the integration of the 
next shell. By dividing the integation into thin shells and taking the limit as the shell width 
Ak goes to zero, eq. (^) can be made self-consistent allowing a non-perturbative solution of 
the path integral. This procedure is equivalent to replacing the bare potential, V, and bare 
wavefunction renormalization, Z, with their k dependent counterparts Up^ and Zp^ on the 
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FIG. 2. Effective critical exponent 7 e g- as function of (a) log 10 (T — T c (k )) and (b) log 10 ((T— 
T c {k )) /k],^ 3 ) at <3? 8 level of truncation. The different symbols correspond to different values of k Q . 
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u /3,k 



In 



p|=* 



sinh 



Z 



(i) 
d.k 



' P,kJ 



2(3 



E 




P 2 /3 + u n + c/: 



+ W 



(6) 



with w p $ now given by u Pi $ = p 2 + L r ^ fc ($)) 1 / 2 . These equations are exact to all loop 
orders since higher-loop contributions are suppressed by additional powers of A A;, allowing 
us to sum over all non-overlapping graphs including daisy and super-daisy diagrams 

Since eq. fl6|) contains two scales, k and (3, there are two ways to proceed at this stage. 
The first would be to fix k and differentiate @ with respect to the temperature; we call this 
the temperature RG 0. We will not discuss this method here but instead will concentrate 
on the other scale, k, which produces the momentum RG. After performing the Matsubara 
sums ]7J in (Q) and differentiating the coupled integral equations with respect to the k we 
have 

dU Ptk _ S d k d \ ( pu Pi $ 



k 



dk 
dk 



sinh 



2 S d k a 



4Sul 9 



-Zk 2 + 9U, 



P,k) 



-n btk (l + (3u kj $(l + n btk )) 



+ 



n b)k (l + n b>k )p 2 u 2 k 



-Zk 2 + 3C4' fe )(l + 2n btk ) - ^Zk 2 (3u k ^(l + 6n btk + 6n b , k 2 ) } 



where n b , k = (e^*.* - l)" 1 and S d = 2/(Air) d / 2 T(d/2). With the initial condition E//j,a($) = 
V($) these equations allow us to integrate down from k = A to k = for a fixed external 
temperature keeping track of not only the relevant but the irrelevant coupling constants as 
well. 



, (7) 
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0^ r ' 1 ^ ~ k am 




(2) 


1.998 ± 0.005 = Tefr/^eff 


4 


0.997 ± 0.004 = Ccff/^cff 


(6) 


0.002 ± 0.004 


(8) 


-1.001 ± 0.005 


(10) 


-2.002 ± 0.005 


<$> 


0.51 ± 0.01 = /Waff 



TABLE I. k dependence of vertex functions at the critical point neglecting the effects of wave- 
function renormalization. 



Typically, the first step made in solving equations like (0) is to expand Up tk as 

oo J 2m ) 

= z #^> 2m , & m) = ^ 2 r } (o), (8) 

and solve the coupled ordinary differential equations for the coefficients gp 2 ™^ after a trun- 
cation of the series. However, it is also possible to solve the coupled equations (|^) without 
making any polynomial expansion or truncation. To accomplish this we simply discretize 
Up t k and Zp t k in $ and k and use finite differences to advance the blocked potential from 
the bare form to the renormalized form. It is desirable to solve ([?]) without polynomial 
truncation since in general the potential will be a non-analytic function of $ at the critical 
point as well as to provide a consistency check of the polynomial expansion. 

In Figure [I] we have plotted the temperature dependence of the renormalized mass, 
coupling constant, and expectation value, $, for d — 3. The Figure shows that all three 
of these quantities vanish at T c continuously indicating a second-order phase transition. 
Exactly how these quantities vanish determines the critical exponents. Since we have control 
over both the temperature and k, we can measure the critical exponents in two ways. The 
first method is to integrate eq. (|7|) down to k = k Q and measure the temperature dependence 
of the thermal parameters near T c {k ) which is defined by the condition fJ^ c ^ ko = &p c (k ),k = 
0. Then using the definitions 

X" 1 = 4 >ko ~ |T - T c (k Q )\™ , (3 -> (3 c (k ) , 
_ dZp cyko ($p c<ko ) R _ R(h \ 

Veft — -g^—^ , P — Pc{Ko) , 

$p jko ~\T-T c (k )f°«, P^PciK), 

®Wo),K ~ h 1 ^ , h^0, (9) 

we can determine the critical exponents. The exponents measured in this way will be 
necessarily k dependent, with the true <i-dimensional critical exponents given in the limit 
k Q — > 0. For k Q ^ the critical exponents measured interpolate between the three- and 
four- dimensional values. 

Another way to determine the critical exponents is to fix T at the true critical temper- 
ature T c = T c (k = 0) and use the k Q dependence of the parameters and finite-size scaling 
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FIG. 3. uj^l as a function of <1> close to T c (T — T c (k Q ) ~ 10 -5 ) obtained by solving the RG 
equations without polynomial truncation. 



laws. Finite size scaling states that near T c the scaling variable is Y = where £ is the 
correlation length of the system and I is the size of the finite dimension. In our analysis there 
are two finite length scales, the blocking length, k' 1 , and the inverse critical temperature, 
(3 C . The scaling variables associated with these two scales are Yp = \T — T c (k )\/k^ UcS and 
Yk = {f3 c k ) l / UcS = /3l^ Ucft , respectively, so that in addition to the standard scaling variable Yg 
we get finite size scaling also in the imaginary time direction governed by the dimensionless 
quantity j3 c . 

Using finite-size scaling we fix T = T c in eq. ([]) and replace \T C — T c {k )\ with k^ UcB 
giving 

= ~ , ^PcJko ~ • (10) 

After extracting u c s by measuring the k a dependence of T c (k ) (limfc o ^o \T c (k )—T c \ ~ kl' Vett ) 
we can determine the effective critical exponents jeg and /3 e s from finite size data. The 
results obtained for the finite size scaling exponents are listed in Table |. Together with a 
measurement of u e s these exponents enable us to calculate 7 e ff, P e g, and Coff- 
in Table || we list the critical exponents we calculated along with previous calculations 
and experimental values. The values obtained using a polynomial truncation of the blocked 
potential agree with those obtained in previous works []5],|6| ; however, the poor agreement with 
experiment suggests that a polynomial truncation of the blocked potential is not appropriate. 
This problem has been pointed out by Morris as well pjfl. We also observe that the (2m)- 
point functions with m > 3 diverge with alternating signs ~ ( — l) (m + 1 )A; m ~ 3 ( 1 +'? e ff). 

This pattern of divergence has been reported elsewhere 0] and provides another indication 
of a non-polynomial structure for the critical effective potential. 

As further evidence for this non-analyticity we can examine the critical effective potential 
obtained from our non-truncated code. In Fig. [3] we have plotted the fourth derivative of 
the blocked potential, uj^J,, near T c for various values of k . From this Figure we see that 

uj^ iiH0 is not identically zero for all $ as would be predicted by a fourth order polynomial 
truncation, i.e. Up ct k = /i^ c fc $ 2 /2 + A^ ci fc$ 4 /4!. Instead, higher order polynomial terms must 

be introduced in order to describe the shape of U$ k _Q. For example, a sixth order term 
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TABLE II. Critical exponents as function of the level of truncation along with the best calcu- 
lations to date and experimental values. In the four-dimensional limit mean field exponents are 
obtained at all levels of truncation. NT indicates results obtained without polynomial trunca- 



tion of the blocked potential. Calculations and experimental values are taken from |TJ and [1C], 
respectively. 



Measured Critical Exponents for d e s = 3 
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1.054 ± 0.002 


0.527 ± 0.003 
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1.171 ± 0.004 


0.585 ± 0.003 
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1.345 ± 0.004 


0.672 ± 0.004 
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1.504 ± 0.004 


0.744 ± 0.004 








NT 


1.234 ± 0.01 




0.315 ± 0.007 


0.0000 


4.65 ± 0.1 


NT + Z 


1.241 ± 0.008 




0.321 ± 0.008 


0.036 ± 0.005 


4.63 ± 0.1 


Best Calculations to Date 


e 5 


1.2390 ± 0.0025 


0.6310 ± 0.0015 


0.3270 ± 0.0025 


0.0375 ± 0.0025 


4.814 ± 0.015 


Experimental Values 




1.23-1.25 


0.624-0.626 


0.316-0.327 


0.016-0.06 


4.6-4.9 



would give Ug c \_Q a parabolic shape near the origin, and higher order terms with increasing 
magnitude and oscillating signs are needed to fit the plateau at large $. 

We have demonstrated that it is possible to determine the critical exponents of the 
scalar A0 4 theory through numerical integration of RG flow equations without polynomial 
expansion. This method provides a straightforward alternative to the e expansion and can 
be used to study field theories in arbitrary dimensions. We have shown through finite-size 
scaling arguments that the effective critical exponents interpolate smoothly between their 
(d + 1)— and d— dimensional values. Our results indicate that the effective potential is a 
non- analytic function of the background field at the critical point, casting serious doubt on 
the validity of a polynomial expansion of the blocked potential in the vicinity of the phase 
transition. In future work we will calculate O(N) exponents and perform studies in two 
dimensions in order to compare our numerical solutions with exact analytic solutions. 
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